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A PROBABILISTIC APPROACH TO DIRICHLET PROBLEMS OF 
SEMILINEAR ELLIPTIC PDES WITH SINGULAR COEFFICIENTS 

By Tusheng Zhang 

University of Manchester 

In this paper, we prove that there exists a unique solution to the 
Dirichlet boundary value problem for a general class of semilinear sec- 
ond order elliptic partial differential equations. Our approach is prob- 
abilistic. The theory of Dirichlet processes and backward stochastic 
differential equations play a crucial role. 

1. Introduction. In this paper, we will use probabilistic methods to solve 
the Dirichlet boundary value problem for the semilinear second order elliptic 
PDE of the following form: 

Au(x) = —f(x,u(x),Vu(x)), Vx G D, 
u(x) \qd = <p, Vx € dD, 

where D is a bounded domain in R d . The operator A is given by 

(1.2) An =1^^ L 3 {x) £\ + £ &,(*) ^ - "div(H" + q(x)u, 

i,j=l 1 ^ i = \ 1 

where a = (aj j(x))i<i j<d : D — > R dxd [d > 2) is a measurable, symmetric ma- 
trix-valued function satisfying a uniform elliptic condition, b = (&i, &2j • • • > bd), 
b = (b±, 02, ■ ■ ■ , bd) ■ D — > R d and q:D^R are merely measurable functions 
belonging to some L p spaces, and /(•, •, •) is a nonlinear function. The oper- 
ator A is rigorously determined by the following quadratic form: 

Q[u,v) = {-Au,v) L 2 {Rd) 

1 f du dv f du 

(L3) =2^J R M x) ^^ dx -^J R M x) ^ {x)dx 



(1.1) 
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(9i> f 
bi(x)u— — dx — / q(x)u(x)v(x) dx. 
D OXi J D 



We refer readers to [14, 18] and [24] for details of the operator A. 

Probabilistic approaches to boundary value problems of second order dif- 
ferential operators have been adopted by many people. The earlier work went 
back as early as 1944 in [15]. See the books [1, 7] and references therein. If 
/ = (i.e., the linear case), and moreover b = 0, the solution u to problem 
(1.1) can be solved by a Feynman-Kac formula 



u(x) = E. 



exp^ T ° q(X(s))ds^(X(r D )) 



for x £ D, 



where X(t), t > is the diffusion process associated with the infinitesimal 
generator 



i,j=l J i=l 



_d_ 

dxi 



to is the first exit time of the diffusion process X(t),t > from the domain 
D. Very general results are obtained in the paper [6] for this case. When b ^ 
0, "div(6-)" in (1.2) is just a formal writing because the divergence does not 
really exist for the merely measurable vector field b. It should be interpreted 
in the distributional sense. It is exactly due to the nondifferentiability of b, all 
the previous known probabilistic methods in solving the elliptic boundary 
value problems such as those in [1, 6, 15] and [13] could not be applied. 
We stress that the lower order term div(fc-) cannot be handled by Girsanov 
transform or Feynman-Kac transform either. In a recent work [5], we show 
that the term b in fact can be tackled by the time-reversal of Girsanov 
transform from the first exit time td from D by the symmetric diffusion X° 
associated with Lq = | Yli j=i ^~( a ij( x ) a§")' symmetric part of A. The 
solution to equation (1.1) (when / = ) is given by 



u(x) = El 



^pT°(T D ))exp / ((a'^X^s^dM^s)) 



o 



(1.5) 



TD 



+ ((a-'b)(X»( S )),dM»(s)))or TD 
1 f TD 
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o 



(b- S)a _1 (6 -b)*(X°(s))ds 

+ r q(X°(s))ds 
Jo 
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where M°(s) is the martingale part of the diffusion X°, r% denotes the reverse 
operator, and (•, •} stands for the inner product in R d . 

Nonlinear elliptic PDEs [i.e., / 7^ in (1.1)] are generally very hard to 
solve. One can not expect explicit expressions for the solutions. However, in 
recent years backward stochastic differential equations (BSDEs) have been 
used effectively to solve certain nonlinear PDEs. The general approach is 
to represent the solution of the nonlinear equation (1.1) as the solution of 
certain BSDEs associated with the diffusion process generated by the linear 
operator A. But so far, only the cases where b = and b being bounded 
were considered. The main difficulty for treating the general operator A in 
(1.2) with b ^ 0, q 7^ is that there are no associated diffusion processes 
anymore. The mentioned methods used so far in the literature ceased to 
work. Our approach is to transform the problem (1.1) to a similar problem 
for which the operator A does not have the "bad" term b. See below for 
detailed description. 

There exist many papers on BSDEs and their applications to nonlinear 
PDEs. We mention some related earlier results. The first result on probabilis- 
tic interpretation for solutions of semilinear parabolic PDE's was obtained 
by Peng in [19] and subsequently in [21]. In [8], Darling and Pardoux ob- 
tained a viscosity solution to the Dirichlet problem for a class of semilinear 
elliptic PDEs (through BSDEs with random terminal time) for which the 
linear operator A is of the form 



where a™ 6 C?(D) and b € C^(D). BSDEs associated with Dirichlet pro- 
cesses and weak solutions of semi-linear parabolic PDEs were considered by 
Lejay in [16] where the linear operator A is assumed to be 



for bounded coefficients a and b. BSDEs associated with symmetric Markov 
processes and weak solutions of semi-linear parabolic PDEs were studied by 
Bally, Pardoux and Stoica in [2] where the linear operator A is assumed 
to be symmetric with respect to some measure m. BSDEs and solutions of 
semi-linear parabolic PDEs were also considered by Rozkosz in [23] for the 
linear operator A of the form 
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Now we describe the contents of this paper in more details. Our strategy 
is to transform the problem (1.1) by a kind of /i-transform to a problem 
of a similar kind, but with an operator A that does not have the "bad" 
term b. The first step will be to solve (1.1) assuming 6 = 0. In Section 2, 
we introduce the Feller diffusion process (n,T,J r t ,X(t),P x ,xeR d ) whose 
infinitesimal generator is given by 

i,j=l J i=l 

In general, X(t), t > is not a semimartingale. But it has a nice martingale 
part M(t), t>0. In this section, we prove a martingale representation the- 
orem for the martingale part M(t), which is crucial for the study of BSDEs 
in subsequent sections. In Section 3, we solve a class of BSDEs associated 
with the martingale part M(t), i > 0: 

(1.7) Y(t) = C+ [ T f(s,Y(s),Z(s))ds- [ T (Z(s),dM(s)). 



The random coefficient f(t,y,z,uj) satisfies a certain monotonicity condition 
which is particularly fulfilled in the situation we are interested. The BSDEs 
with deterministic terminal time were solved first and then the BSDEs with 
random terminal time were studied. In Section 4, we consider the Dirichelt 
problem for the second order differential operator 

i,j=l J i=l 

where b{ € LP for some p> d and q € L@ for some /3 > | . We first solve the 
linear problem with a given function F 

L2U = F, in D, 
u = tp, on dD, 



(1.9) 

and then the nonlinear problem 
(1.10) 



L2U = — g(x, u(x), Vn(x)), in D, 
u = (p, on dD, 

with the help of BSDEs. Finally, in Section 5, we study the Dirichlet problem 

J Au(x) = — f(x,u(x)), VxGD, 

\u(x)\ dD = ip, yxedD, 

where A is a general second order differential operator given in (1.2). We 
apply a transform we introduced in [5] to transform the above problem to 
a problem like (1.10) and then a reverse transformation will solve the final 
problem. 



(1.11) 
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2. Preliminaries. Let A be an elliptic operator of the following general 
form: 

A = \*t +E¥x)|:-»div(5r + 

i,j=l J i=l 

where a = {a^x)) : D -> i? dxa! (d > 2) is a measurable, symmetric matrix- 
valued function which satisfies the uniform elliptic condition 

d 

(2.1) A|£| 2 < < A|£| 2 , G i? d and x £ D 

for some constant A, A > 0, b = (b\, . . . ,b<i),b = (b\, . . . ,bd) D ^ R d and : 
D R are measurable functions which could be singular and such that 

|6| 2 GL P (D), |6| 2 GLf(D) and q€L p (D), 

for some p > |. Here D is a bounded domain in whose boundary is 
regular, that is, for every x G dD, P(t^ = 0) = 1, where r|j, is the first 
exit time of a standard Brownian motion started at x from the domain D. 
Let / : R d x R x — > R be a measurable nonlinear function. Consider the 
following nonlinear Dirichlet boundary value problem: 

, s ( Au(x) = -f(x,u(x), Vu(x)), VxG-D, 

K ' \u(x)\ dD = (p, VxedD. 

Let W 1,2 {D) denote the usual Sobolev space of order one: 
W l ' 2 {D) = {u G L 2 {D) : Vu G L 2 {D; R d )}. 

Definition 2.1. We say that u G W 1,2 (D) is a continuous, weak solu- 
tion of (2.2) if: 

(i) for any 0G W 1,2 (D), 

(j) dx 

2 ' In lJK ' dxi dxi ^—f In ' dxj 



y~] [ bi{x)u^-dx- [ q(x)u(x)<p dx = [ f(x, 
~^JD OXi J D Jn 



u, Vu)4>dx, 



(ii) ueC(D), 

(hi) lim y -> x u{y) = tp(x), Vx G dD. 

Next we introduce two diffusion processes which will be used later. 
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Let {Q, J", T t ,X(t),P x ,x € R d ) be the Feller diffusion process whose in- 
finitesimal generator is given by 




where Tt is the completed, minimal admissible filtration generated by X(s), 
s > 0. The associated nonsymmetric, semi-Dirichlet form with L\ is defined 
by 

Qi(u,v) = (-Liu,v) L 2 iRd) 

(2-4) 

= -V/ a ij (x) — — dx-S2 bi{x)—v{x)dx. 

The process X(t), t > is not a semimartingale in general. However, it 
is known (see, e.g., [6, 10, 12] and [17]) that the following Fukushima's 
decomposition holds: 

(2.5) X(t) = x + M(t) + N(t) Pa.-a.s-, 

where M(t) is a continuous square integrable martingale with sharp bracket 
being given by 

(2.6) {{M\M%= f ai!j (X(s))ds, 

Jo 

and N(t) is a continuous process of zero quadratic variation. Later we also 
write X x (t), M x (t) to emphasize the dependence on the initial value x. Let 
M denote the space of square integrable martingales w.r.t. the filtration 
Tti t>0. The following result is a martingale representation theorem whose 
proof is a modification of that of Theorem A. 3. 20 in [12]. It will play an 
important role in our study of the backward stochastic differential equations 
associated with the martingale part M. 

Theorem 2.1. For any L € M., there exist predictable processes Hi(t),i = 
1, . . . , d such that 

d ft 

(2.7) L t = Y J Hi^dM^s). 

i=i Jo 

Proof. It is sufficient to prove (2.7) for < t < T, where T is an arbi- 
trary, but fixed constant T. Recall that JVl is a Hilbert space w.r.t. the in- 
ner product {K\ , K2)m = E[{{K\ , K2))t] , where ((K\ , K2}} denotes the sharp 
bracket of K\ and K2. Let J\4.\ denote the subspace of square integrable mar- 
tingales of the form (2.7). Let R a ,a > be the resolvent operators of the 
diffusion process X(t), t>0. Fix any g € Cb(R d ), we know that R a g € D(L±) 
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and L\R a g = ctR a g — g. Moreover, it follows from [12] and [17] that 

R a g(X(t))-R a g(X(0)) = f (V R a g(X(s)),dM(s)) 

Jo 

+ f (aR a g-g)(X(s))ds. 
Jo 

Hence, 

Jf-= f e- as (VR a g(X(s)),dM( S )) 
Jo 

= e- at R a g(X(t)) - R a g(X(0)) + f e - as g(X(s)) ds 

Jo 

is a bounded martingale that belongs to M\. The theorem will be proved 
if we can show that = {0}. Take K G Mi - Since Mi is stable under 
stopping, by Lemma 2 in Chapter IV in [22], we deduce ((K,L)) =0 for all 
L € M.\. In particular, ((K, J)) = 0. From here, we can follow the same proof 
of Theorem A.3.20 in [12] to conclude K = 0. □ 

We will denote by (fi, T , J=f, X°(t),P®,x € R d ) the diffusion process gen- 
erated by 




The corresponding Fukushima's decomposition is written as 

X°(t) = x + M°(t) + N°(t), t>0. 

For v G W 1)2 {R d ), the Fukushima's decomposition for the Dirichlet process 
v(X°(t)) reads as 

(2.9) v(X°(t)) = v(X°(0)) + M v (t) + N v (t), 

where M v (t) = £ Vv(X° (a)) • dM° (s) , N v (t) is a continuous process of zero 
energy (the zero energy part). See [3, 4, 12] for details of symmetric Markov 
processes. 

3. Backward SDEs with singular coefficients. Let (CI, J 7 , J^) be the prob- 
ability space carrying the diffusion process X(t) described in Section 2. Re- 
call M(t), t > is the martingale part of X. In this section, we will study 
backward stochastic differential equations (BSDEs) with singular coefficients 
associated with the martingale part M(t). 

3.1. BSDEs with deterministic terminal times. Let f(s,y,z, co) : [0,T] x 
R x R d x Cl — > R be a given progressively measurable function. For simplicity, 
we omit the random parameter to. Assume that / is continuous in y and 
satisfies: 
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(A.l) (y 1 -y 2 )(f{s,y 1 ,z) - f{s,y 2 ,z)) < - y 2 \ 2 , 

(A.2) \f(s,y,zi) - f{s,y,z 2 )\ < d 2 \z 1 - z 2 \, 



(A.3) \f(s,y,z)\ < \f(s,0,z)\ + K(s)(l + \y\), 



where di(-), K(s) are a progressively measurable stochastic process and d 2 
is a constant. Let £ € L 2 (0, .Ft, -P)- Let A be the constant defined in (2.1). 



Theorem 3.1. Assume E[e~^ 2di(a)<fe|£|2] < ^ E[Jq K (s) ds] < oo and 



Then, i/iere exists a unique (Tt- adapted) solution (Y,Z) to the following 



BSDE: 

(3.1) Y(t) = £ + [ T f(s, Y(s), Z(s)) ds - [ T (Z(s),dM(s)), 
Jt Jt 

where Z(s) = (2i(s), . . . , Z d {s)). 

Proof. We first prove the uniqueness. Set d(s) = — 2d\(s). Suppose 
(Y 1 ^),^ 1 ^)) and (Y 2 (t) , Z 2 (t)) are two solutions to equation (3.1). Then 

d(\Y\t)-Y 2 (t)\ 2 ) 

= -2(Y\t)-Y 2 (t))(f(t,Y 1 (t),Z\t))-f(t,Y 2 (t),Z 2 (t)))dt 

(3.2) 

+ 2(Y\t) - Y 2 (t)){Z\t) - Z 2 (t),dM(t)) 
+ (a(X(t))(Z l (t) - Z 2 {t)),Z\t) - Z 2 (t)) dt. 



By the chain rule, using the assumptions (A.l), (A.2) and Young's inequality, 




we get 



|Y x (t) - Y 2 (t)\ 2 eti d ^ ds 




x (f(s, Y\s), Z\s)) - f(s, Y 2 (s),Z 2 (s))) ds 
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(3.3) < - jfVo <*(«)<*« |yi( s ) -Y 2 (s)\ 2 d(s)ds 

-2 £ efo ^ | Y 1 (s) -Y 2 {s)\ 2 d 1 (s) ds 

+ 2 j\^ di > u)du d 2 \Y\s)-Y 2 {s)\\Z 1 {s)- Z 2 {s)\ds 

-2 j\fo d ^ du (Y\s)-Y 2 (s)){Z 1 (s)-Z 2 (s),dM(s)) 

rT 

<C\J^ e^ d ^ du \Y 1 {s)-Y 2 {s)\ 2 ds 

+ \ [ eK d ^(a{X{s)){Z\s) - Z 2 (s)), (Z\s) - Z 2 {s))) ds 

-2 j\fo d ^ du (Y 1 (s)-Y 2 (s))(Z 1 (s)-Z 2 (s),dM(s)). 

Take expectation in above inequality to get 

E[\Y\t) - Y 2 {t)\ 2 eti d ( s)ds ] < C X J T E[ef° d ( u ) du \Y\s) - Y 2 (s)\ 2 ] ds. 

By Gronwall's inequality, we conclude Y l {t) = Y 2 {t) and hence Z l {t) = 
Z 2 (t) by (3.3). 

Next, we prove the existence. Take an even, nonnegative function <f> € 
Cfi°(R) with f R <j>(x) dx = l. Define 

fn(t,y,z)= f(t,x,z)4> n (y-x)dx, 
JR 

where 4> n (x) = ruj)(nx). Since / is continuous in y, it follows that f n (t, y, z) — > 
f(t,y,z) as n — > oo. Furthermore, it is easy to see that for every n > 1, 

(3.4) \fn(t,yi,z) - f n (t,y 2 ,z)\<C n \y 1 -y 2 \, yi,y2^R, 
for some constant C n . Consider the following BSDE: 

(3.5) Y n (t)=C+ [ T f n (s,Y n (s),Z n (s))ds- [ T (Z n (s),dM(s)). 

Jt Jt 

In view of (3.4) and the assumptions (A. 2), (A. 3), it is known (e.g., [20]) 
that the above equation admits a unique solution (Y n ,Z n ). Our aim now is 
to show that there exists a convergent subsequence (Y nk ,Z nk ). To this end, 
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we need some estimates. Applying Ito's formula, in view of assumptions 
(A.1)-(A.3) it follows that 

\Y n {t)\ 2 eJ» d ^ ds + [ T ef° d ^ du (a(X(s))Z n (s),Z n (s))ds 



+ 2 ^ efo rf H du Y n (s)f n (s, Y n {s),Z n {s)) ds 



T 



>J° d M du Y n (s)(Z n (s),dM(s)) 



<\£f e fod{s)ds _ / e So d ^) du Y 2 {s)d{s)ds 



T 



T 



^ d{u)du d l (s)Y 2 (s)ds + 2C / e^ d ^ u)du \Y n {s)\\Z n {s)\ds 



+ 2 / e f° d ^ du \Y n (s)\f(s,0,0)ds 



(3.6) 



T 



^ d ^ du Y n (s)(Z n (s),dM(s)) 



T 



< |f |2 e / d{s)ds + Cx / g / s d{u)du Y 2^ dg 



e% d M du (a{X{s))Z n {s),Z n {s))ds 

+ J e^ d ( u ) du Y 2 {s)ds + J e^ d ( u ) du \f( s ,0,0)\ 2 ds 

^ d{u)du Y n (s){Z n (s),dM(s)). 
Take expectation in (3.6) to obtain 

E[\Y n (t)\ 2 eti d ^ ds ] + ^E j\fo d M d -(a(X(s))Z n (s),Z n (s))d S 

(3.7) <E[\Z\ 2 eSo d ( s ) ds } + C x J T E[efo d ( u ) du Y*(s)}ds 



+ E 



T 



>J° d M du \f(s,0,ty\ 2 ds 
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Gronwall's inequality yields 



(3.8) 



sup sup E[\Y n (t)\ 2 efo 

n 0<t<T 



<C^E[\^\ 2 e^ d{s)ds ]+E j\-fo d ^ du \f( s ,0,0)\ 2 ds 



and also 



(3.9) 



supE 



Jo d M d -( a (X(s))Z n ( S ),Z n (s))ds 



< oo. 



Moreover, (3.6)-(3.9) further imply that there exists some constant C such that 
E 



sup Y 2 {t)e& d ^ ds 

0<t<T 



<C + CE 

<C + CE 
<C + CE 



sup 

0<t<T JO 



^ d{u)du Y n (s)(Z n (s),dM(s)) 



7 \ 2 ^ d{u)du Y 2 { S )(a{X{s))Z n {s),Z n {s))d S X 



(3.10) 



<C + -E 



+ C X E 

In view of (3.9), this yields 
(3.11) sup£ 



o 



sup {eWtt d ^ du \Y n (s)\) 

0<s<T 



1 e fo d M d u(a(X(s))Z n (s),Z n (s))ds^ 

o 



sup (e^W^(s)) 

0<s<T 



fo d ^ du (a(X(s))Z n (s),Z n (s))ds 



sup Y 2 {t)e^ d ^ s)ds 

0<t<T 



< OO. 



By (3.9) and (3.11), we can extract a subsequence such that Y nfc (i)e( 1//2 ) d ( s ) ds 
converges to some Y(t) in L 2 (Q,, L°°[0,T]) equipped with the weak star 
topology and Z nk [t)e^/ 2 ^o d i s ) ds converges weakly to some Z(t) in L 2 (VLt] R), 
where = [0,T] x Q. Observe that 

Y nk {t)e^ti d ^ ds 



T 



(i/2)/ ^ (s )^ + / e (imio d (^f nk{SiYnh (s),Z nk (s))ds 



12 
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e WJo d ( u ) du Y nk (s)d(s)ds 
^mSo d ^) du {Z nk {s),dM{s)). 



Letting k — > oo in (3.12), using the monotonicity of /, following the same 
arguments as that in the proof of Proposition 2.3 in Darling and Pardoux 
in [8], we can show that the limit (Y, Z) satisfies 

Y(t) =e (V2)/o T <*W<k£ 
(3.13) 



+ 



,(1/2) J ( ; d(u) du ^ e -(l/2) J S d(u) duf( a y e -{l/2) / 3 d(u) **Z{a)) ds 



T 



Y(s)d(s)ds 



(Z(s),dM(s)). 



Set 



Y(t) = e" (1/2) Jo <*(") du y (t), Z(t) = e~ {1/2) Jo *(«) . 



An application of Ito's formula yields that 
Y(t) = €+ C f(s,Y(s),Z(s))ds 



(Z(s),dM(s)), 



namely, (Y,Z) is a solution to the backward equation (3.1). The proof is 
complete. □ 

3.2. BSDEs with random terminal times. Let f(t,y,z) satisfy (A.l)- 
(A.3) in Section 3.1. In this subsection, set d(s) = —2d\{s) + bd\. The fol- 
lowing result provides existence and uniqueness for BSDEs with random 
terminal time. Let r be a stopping time. Suppose £ is J-~ r -measurable. 

Theorem 3.2. Assume E[e$o d{ ~ s)ds \£\ 2 ] < oo, E[JqK(s)cIs] < oo and 



(3.14) 



E 



Jo<W«|/( S)0 ,0)| 2 ds 



< oo, 



for some 5 > j, where X is the constant appeared in (2.1). Then, there exists 
a unique solution (Y, Z) to the BSDE 

(3.15) Y{t) = i+ f f(s,Y(s),Z(s))ds- [ T (Z(s),dM(s)). 
Furthermore, the solution (Y, Z) satisfies 



(3.16) E 



i$o d ^ du Y 2 (s)ds 



< oo, 



E 



Jo d H du \Z(s)\ 2 ds 



< oo, 
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and 

(3.17) E \ sup {e^ d{u)du Y 2 (s)} 

K)<a<T 

Proof. After the preparation of Theorem 3.1, the proof of this theorem 
is similar to that of Theorem 3.4 in [8], where di(s), c?2 were both assumed 
to be constants. We only give a sketch highlighting the differences. For every 
n > 1, from Theorem 3.1 we know that the following BSDE has a unique 
solution (Y n , Z n ) on < t < n: 

rrAn rrAn 

(3.18) Y n (t)=E[£\F n ]+ / f(s,Y n (s),Z n (s))ds- (Z n (s) , dM (s)) . 

JtM JtM 

Extend the definition of (Y n , Z n ) to all t > by setting 

Y n {t) = E[£\JF n ], Z n {t) = for t > n. 

Then the extended (Y n ,Z n ) satisfies a bsde similar to (3.18) with / replaced 
by A{s<nAr}/( s ) Di z )- Let n>m. By Ito's formula, we have 

\Y n (tAT)-Y m (tAT)\ 2 eti ATd ^ ds 
pnAT 

+ / efo d ^ d -(a(X(s))(Z n (s) - Z m (s) X{s<mAT} ), 

JtAT 

Z n (s) - Z m (s)x{ s<mA 

r})ds 

(3.19) =e^ T ^)^{E[i\F n ]-E[i\F m ]f 
e C T d(u)du\y n ( s /\ r ) _ Y m (s A T)\ 2 d(s) ds 



tAT 

rnAr 

+ 2/ efo Td ^ du (Y n (sAT)-Y m (sAr)) 

JtAT 

x (f(s,Y n (sAT),Z n (sAT)) 

- /(s, Y m (s A r), Z m (s A r)))ds 

e J^ d(u)du {Yn{s A r) _ Ym{s A r)) 

mAr 

x f{s,Y m (s Ar),Z m {s Ar))ds 



e 

tAT 

mAr 



s At 

tAT 



f° ATd M du (Y n (s At)- Y m (s A r)){Z n {s A r),dM(s)) 

j ( rd(u)du {Yn{s A r) _ Ym{s A T)){Zm{s A T ) )dM ( s )). 



14 T. ZHANG 

Choose Si, 62 such that j < 5\ < 5 and < 82 < 5 — 5±. In view of the (A.l) 
and (A. 2), we have 

fflAT 

2/ e/« Td W 4 (7 tt ( S Ar)-7 m ( S Ar)) 

x (/(a, y n (s A r), Z n (s A r)) - /(a, F m (s A r), Z m (s A r))) (is 



<_ 2 / e /o Td ^ du (Y n (sAT)-Y m (sAr)) 2 d 1 (s)ds 
(3.20) ^ 

f nAr A 

+ / e-To* Td W^(y n ( S Ar)-^ n ( S Ar)) 2 da 



tAr 
nAr 



1 rnivr 

+ TT / e /Wd >W*))(^«(s) " ^™(s)X {s <™At}), 



JWr 

Z n (s) - Z m (s)x{ s <mAr}) 

On the other hand, by (A. 3), it follows that 



e / s 'dMdu(Y n ( a a r) — F m (a A r))/(a, Y m (s A r), Z m (a A r)) ds 

mAr 

/•nAr 

(3.21) <<5 2 d|/ e-K Td ^ du (Y n (sAT)-Y m (sAT)) 2 ds 

e Io AT ^)^(\f( s ,0,0)\+K + K\E[^ m ]\) 2 ds. 



tAr 



1 

Take expectation and utilize (3.19)-(3.21) to obtain 
E[\Y n (t At) - y,„(t A t)|V"' "<•>"■] 



nAr 

gj 

Mr 



/o<«)^(a(X(a))(Z n (a)-Z m (a)x {s < mA r}), 

Z„(s) - Z m (s)x{ s < m .Ar}) ^ 
f°d(u)du( Yn ( s A r ) _ A r ^2 ds 

^ d{s)ds{mFn] _ mFm])2] 

e C Td M d "(\f( sA 0)\+K + K\E[ti\T m }\) 2 d S 



(3.22) 

r mAr 

+ {5-5 l - 5 2 )d 2 2 E 



tAr 



mAr 



Since the right-hand side tends to zero as n, m — )■ 00, we deduce that 

{( e ( 1 /2)/o tAT ^)*y n (t) |e (V2)/ tA M-)* Zn ( < )) } 
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converges to some (Y,Z) in M 2 (0,t;R x R d ). Furthermore, for every t > 0, 
eC 1 / 2 )^ d ^ ds Y n (t) converges in L 2 . We may as well assume 

(3.23) Y(t)= lim e {l/2) C Td ^) ds Y n {t) 

n—^oo 

for all t. Observe that for any n > t > 0, 

i-nAr 

= e (l/2)/ T d(»)d»£[£prj + / e (l/2)/ S ^)^/( s ,y n ( s ),Z n ( S ))ds 



(3.24) 



rAt 



e (1/2) ^o ATd(n)dM y n (s)d(s)ds 

rAt 
nAr 



'rAt 

Letting n — )• oo yields that 



rAt 



(3.25) g-ti/^/^dW*^^^ 
1 



Y(s)d(s)ds- (Z(s),dM(s)). 

* JrAt JrAt 

Put 

Y(t) = e -i 1 l 2 )C T ^ s ) ds Y{t), Z(t) = e - {1/2) ^ Td{s)ds Z{t). 
An application of Ito's formula and (3.25) yield that 

(3.26) Y(t) = £+ f T f(s,Y(s),Z(s))ds- f (Z(s),dM(s)). 

JrAt JrAt 

Hence, (Y, Z) is a solution to the bsde (3.15) proving the existence. To obtain 
the estimates (3.16) and (3.17), we proceed to get an uniform estimate for 
Y n {s) and then pass to the limit. Let 81,62 be chosen as before. Similar to 
the proof of (3.8), by Ito's formula, we have 

ftiAr 

|F n (tAr)|Vo T Ws+ / eS° d M du (a(X(s))Z n (s),Z n (s))ds 

JtAr 
A rnAr 

<\E[£\F n ]\ 2 efo Td ( s ) ds - / eK d ^ du \Y n {s)\ 2 d{s)ds 

JtAr 

rnAr 

_ 2 / e^^^dx^Y^ds 



16 

+ 2 
+ 2 

(3.27) - 2 



n.Ar 



tAT 
nAT 



tAT 
nAr 



tAr 
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Jo^)^d 2 \Y n ( S )\\Z n (s)\ds 

Jo^)^\ Yn (s)\\f(s,0,0)\ds 
^ d(u)du Y n (s)(Z n (s),dM(s)) 



<\E[t\F n ]\ 2 efo nAT <i(* 



uAt 



uAt 



tAT 



e So d ( u ) du 5d 2 2 Y 2 {s)ds 

fo^)du 5id 2 Y 2^ ds 
i rnAT 

+ / e /o^)^( a (X( S ))Z n ( S ),Z n ( S ))d S 

/•nAT i rnAT 

+ / 524ef° d ^ du Y 2 (s)ds + —» eJo d ( u ) du \f( s ,0,0)\ 2 ds 

JtAT °2«2 ^iAr 

^ d W du y n ( s )(z n ( s ),dM( s )). 



uAt 



tAT 



Recalling the choices of d(s), 5\ and 52 , using Burkholder's inequality, we 
obtain from (3.27) that 



E 



sup 

0<t<n 



\Y n {t Ar)| 2 e^o ATd(s)ds 



(3.28) <E[\i\ 2 e^ Td{s)ds ] + E fVo rf(«)^_L|/( a ,0,0)| 2 ds 

JO 02^2 

e 2 ^ d H d «y n 2 ( s )(a(X( S ))Z n ( S ),Z n ( S ))d S 

< Vo" AT d « + £ T e /o d M d « J_ |/(« 0, 0) | 2 ds 

Jo hd A 2 



+ 2CE 



1/2 



2 



sup \Y n (t At)|V» T 

0<t<n 



nAr 



Jo d (-) d -(a(X( S ))Z n ( S ),Z n ( S ))ds 



In view of (3.27), as the proof of (3.9), we can show that 



(3.29) sup£ 

n Uo 



nAT 



Jo d ^ du (a(X( S ))Z n (s),Z n ( S ))ds 



< oo. 
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(3.29) and (3.28) together with our assumptions on / and £ imply 



(3.30) supE 

n 

Applying Fatou lemma, (3.17) follows. □ 



sup \Y n (tAT)\ 2 eti ATd(s)ds 

0<t<n 



< oo. 



3.3. A particular case. Let f(x, y, z) : R d x R x R d — > R be a Borel mea- 
surable function. Assume that / is continuous in y and satisfies: 

(B.l) (yi - y 2 ){f{x,y 1 ,z) - f(x,y 2 ,z)) < -ci(x)\yi - y 2 \ 2 , where c x (x) is a 

measurable function on R d . 
(B.2) \f(x,y,zt) - f(x,y,z 2 )\ < c 2 \z\ - z 2 \. 
(B.3) \f(x,y,z)\ < \f(xAz)\+c 3 (x)(l + |y|). 

Let D be a bounded regular domain. Define 

(3.31) T%=M{t>0:X x (t)(£D}. 

Given g € Cb(R d )- Consider for each x € D the following BSDE: 



Y x (t)=g(X x (rf ) )) + 



f(X x (s),Y x (s),Z x (s))ds 



(3.32) 



tAr? 



(Z x (s),dM x (s)), 



thri 



where M x (s) is the martingale part of X x {s). As a consequence of Theo- 
rem 3.2, we have the following theorem. 



Theorem 3.3. Suppose c 3 G L P (D) for p > |, 



exp 



(-2ci(X(s)) + 5c^)ds 



< oo, 



/or some <5 > j and 



E , 



|/(X( S ),0,0)| 2 d S 



< oo. 



The BSDE (3.32) admits a unique solution (Y x (t), Z x (t)) . Furthermore, 



(3.33) 



sup (0)| < oo. 



4. Semilinear PDEs. As in previous sections, (X(t),P x ) will denote the 
diffusion process defined in (2.5). 

4.1. Linear case. Consider the second order differential operator 
(4.1) L 2 = \j:^ (a,(x) A) + £ Hx) A. + q{x) . 

i,j=l J i=l 
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Let D be a bounded domain with regular boundary (w.r.t. the Laplace 
operator A) and F(x) a measurable function satisfying 

(4.2) \F{x)\<C + C\q{x)\. 

Take ip € C(dD) and consider the Dirichlet boundary value problem 

L 2 u = F, in D, 
u = <p, on dD. 



(4.3) 



Theorem 4.1. Assume (4-2) and that there exists xq^lD such that 



E 



■i-f) 



exp 



l(X(s))ds 



< 00. 



Then there is a unique, continuous weak solution u to the Dirichlet boundary 
value problem (4-3) which is given by 



(4.4) u(x) = E x 

Proof. Write 



<p{X{r%)) + I ° eti q{x{s))ds F{X{t)dt 



u l {x) = E x [i P {X{r x D ))], 



and 



u 2 (x) = E x 



° eti q{x{s ^ ds F{X{t))dt 



We know from Theorem 4.3 in [6] that u\ is the unique, continuous weak 
solution to the problem 



(4.5) 



L2U = 0, in D, 
u = (f, on dD. 

So it is sufficient to show that u 2 is the unique, continuous weak solution to 
the following problem: 



(4.6) 



L 2 u = F, in D, 
u = 0, on dD. 



By Lemma 5.7 in [6] and Proposition 3.16 in [7], we know that u 2 belong to 
Cq{D). Let G p,/3 > denote the resolvent operators of the generator L 2 on 
D with Dirichlet boundary condition, that is, 



Gf)f(x) = E x 



-^e^^ x ^ ds f{X{t))dt 



By the Markov property, it is easy to see that 

P(v s (x)-PGf}U2(x))=PGf ) F(x). 
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Since Gp is strong continuous, it follows that 

lim p{u 2 - pGau 2 ) = F 

in L 2 (D). This shows that u 2 € D(L 2 ) C W l > 2 (D) and L 2 u 2 = F. The proof 
is complete. □ 

4.2. Semilinear case. Let g(x,y,z) :R d x R x R d — > i? be a Borel mea- 
surable function that satisfies: 

(C.l) (yi -y 2 )(g(x,y 1 ,z) -g(x,y 2 ,z)) < -ki(x)\yi — 2/2 1 2 , 
(C.2) -^(a;,y,2r 2 )| < fc 2 |«i - z 2 \, 

(C.3) 1^,2/^)1 <C + C|g(x)|, 

where k\{x) is a measurable function and k 2 ,C are constants. Consider the 
semilinear Dirichlet boundary value problem 

Uj\ ( L 2 u= -g(x,u(x),Vu(x)), in D, 

\ u = <p, on dD, 

where ip G C(dD). 

Theorem 4.2. Assume 

< oo, 



F x 



I D (q(X(s))-2k 1 (X(s))+5k 2 2 )ds 



for some S > j. 

The Dirichlet boundary value problem (4-7) has a unique continuous weak 
solution. 

Proof. Set f(x, y, z) = q(x)y + g(x, y, z). According to Theorem 3.3, for 
every x € D the following BSDE: 

Y x {t) = <j ) {X x (T x D ))+ / ° f{X x {s),Y x {s),Z x (s))ds 

(48) 

3 (Z x (s),dM x (s)), 

X 

D 

admits a unique solution (Y x (t), Z x (t)),t > 0. Put uo(x) = Y x (0) and vq(x) = 
Z x (0). By the strong Markov property of X and the uniqueness of the BSDE 

(4.8) , it is easy to see that 

(4.9) Y x (t)=u (X x (t)), Z x {t) = v (X x {t)), Q<t<T x D . 

Now consider the following problem: 

L 1 u = -f(x,u (x),v (x)), in D, 
u = cp, on dD, 



(4.10) 
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where L\ is defined as in Section 2. By Theorem 4.1, problem (4.10) has a 
unique continuous weak solution u(x). As u € V^ 1,2 (.D), it follows from the 
decomposition of the Dirichlet process u(X(t ArJ)) (see [12, 17]) that 

u{X{tAr x D ))=^{X x {T x D )) + f T ° f(X x (s),u (X(s AT^M^is AT^ds 



{Vu{X(sAT x D )),dM x (s)) 



(4.11) 



tAr* 



^(X x (r^))+ / f(X x (s),Y x (s),Z x (s)))ds 



(Vu(X{sAT?))),dM x {s)). 



thri 



Take conditional expectation both in (4.11) and (4.8) to discover 
Y x {tAr x D )=u{X{tArf ) )) 



E 



y{X x {T x D ))+ / f(X x (s),Y x (s),Z x (s))ds 



thrf 



In particular, let t = to obtain u(x) = uq(x). On the other hand, comparing 
(4.8) with (4.11) and by the uniqueness of decomposition of semimartingales, 
we deduce that 



{Vu{X{sAT x D )),dM x {s)) 



tAr? 



(Z x (s),dM x (s)) 



iAr? 



for all t. By Ito's isometry, we have 
E 



(4.12) 



E 



E 



((Vu(X(s))-Z x (s)) X{s<T * } ,dM x (s)) 



((Vu(X(s)) - vo(X x (s)))x {s<T *},dM x (s)) 
a(X x (s))(Vu(X(s))-v (X x (s))), 

Vu(X(s)) - D (X :r (s)))x{ s <r- } ds 



2-1 



0. 



By Fubini theorem and the uniform ellipticity of the matrix a(x), we deduce 
that 

P°(\Vu - v \ 2 ) = E[\Vu(X(s)) - v (X x (s))\ 2 X{s<Th} ] = 
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a.e. in s with respect to the Lebesgue measure, where P?h(x) = E x [h(X(t)), 
t < Tq] . The strong continuity of the semigroup P® , s > implies that 

(4.13) |Vu - v \ 2 (x) = lim P?(\Vu - v \ 2 ) = 

S— ¥0 

a.e. Returning to problem (4.10), we see that u actually is a weak solution 
to the nonlinear problem: 

Lqu = — f(x, u(x), Vu(x)), in D, 
u = 4>, on dD. 

Suppose u is another solution to the problem (4.14). By the decomposition 
of the Dirichlet process u(X x (s)), we find that (u(X x (s)),Vu(X x (s))) is 
also a solution to the BSDE (4.8). The uniqueness of the BSDE implies 
that u(X x (s)) = Y x (s). In particular, u(x) = uq(x) = Y x (0). This proves the 
uniqueness. □ 



(4.14) 



(5.1) 



5. Semilinear elliptic PDEs with singular coefficients. In this section, we 
study the semilinear second order elliptic PDEs of the following form: 

J Au(x) = —f(x, u(x)), Vx € D, 
\u(x)\ dD = (p, VxedD, 

where the operator A is given by 

i,j = l i=l 

as in Section 2 and <p G C(dD). Consider the following conditions: 

(D.l) ( yi - y 2 )(f(x, yi ) - f(x,y 2 )) < -Ji(x)| yi - y 2 \\ 
(D.2) \f(x,y,z)\<C, 

where J\(x) is a measurable function, C is a constant. The following theorem 
is the main result of this section. 

Theorem 5.1. Suppose that (D.l), (D.2) hold and 



cxp 



((a^ 1 6)(X ( S )),dM°( S )) 

+ ^j T \{a- l b){X Q (s)),dM\s))^or 
1 f' TD 

(5.2) -=- / (6-6)a- 1 (6-6)*(X°(s))ds 



2 



[TD fTD 

+ q(X°(s))ds-2 Ji(X°(s))ds 
Jo Jo 



< oo 
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for some x G D, where X° is the diffusion generated by Lq as in Section 2 
and td is the first exit time of X° from D. Then there exists a unique, 
continuous weak solution to equation (5.1). 

Proof. Set 

Z t = e X p|^V^ 1 &)(^^s)),dM ( S )) + ^((«"^)(^ (s))^M°(^))) or t 
(5.3) - \ I (b - b)a~ l (b - b)*(X°(s)) ds + f q(X°(s)) ds 



2 



*J 1 (Z°( S ))d S [ > . 



Put 



ft 

M(t)= l ((a^ 1 b){X (s)),dM°(s)) fort>0. 



Let R > so that D C B R := B(0,R). By Lemma 3.2 in [5] (see also [9]), 
there exits a bounded function v S Wq ,p (Br) C Wq' 2 {Br) such that 

(M(t))or t = -M(t) + N v (t), 

where iV" is the zero energy part of the Fukushima decomposition for the 
Dirichlet process v(X°(t)). Furthermore, v satisfies the following equation 
in the distributional sense: 

(5.4) div(aVv) = -2div(S) in B R . 

Note that by Sobolev embedding theorem, v G C{R d ) if we extend v = on 
D c . This implies that M and N v are continuous additive functionals of X° 
in the strict sense (see [9, 12]), and so is t — )• (M(t)) o rt- Thus, 

\(a- 1 b)(X i> (8)),dM (s)))or t 



< ((a- 1 6)(X°( S )),dM°( S ))+^(t) 
t 

((a" 1 6)(X°(s)), dM°(s)) + v(X°(t)) - v(X°(0)) - M v (t) 

{(a-^X^s)), dM°(s)) + v(X°(t)) - v(X°(0)) 
t 

(Vv(X {s)),dM°(s)). 



A PROBABILISTIC APPROACH TO DIRICHLET PROBLEMS OF PDES 23 



Hence, 

e «(*°(*)) 
Zt = e v(XO(0)) 



xexp^ (a- 1 (b-b-aVv)(X (s)),dM°(s)) -2 J J X (X° (s)) ds 
+ J (v-^( b -b-aVv)a- 1 (b-b-aVv)*^{X (s))ds 



Note that Zt is well defined under P® for every x € -D. Set /i(x) = e"^. 
Introduce 

i = 5 E Sr. ( a ^ Sr) + ~ to ~ ( a ^Ux)] ^ 

i,j=l % V 3 J i=l % 

- (b- b,Vv){x) + -(Vv)a(Vv)*(x) + q(x). 

Let (Q,J-,J-t,X(t),P x ,x € R d ) be the diffusion process whose infinitesimal 
generator is given by 

L = \ S Sr. \ ai ^ Sr) + E^( x ) - to ~ ( a Vv)i(x)} 

i,j=l 1 V j/ i=l 4 

It is known from [17] that P x is absolutely continuous with respect to P x 
and 



(5.5) 



dP x 



= Z t , 



where 



(5.6) 



Z t = «xp(J ( a ^(b -b- aVv)(X (s)),dM°(s)) 

-J Q(6 - b - aVw)a" 1 (6 - b - aVv)*^j (X°(s)) ds\ . 



Put 

f(x,y) = h(x)f(x,h~\x)y). 

Then 

(yi-y2)(f(x,yi) - f(x,y 2 )) < -Ji(x)\yi -y 2 \ 2 - 
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Consider the following nonlinear elliptic partial differential equation: 
(5.7) 

In view of (5.5), condition (5.2) implies that 



j Au(x) = — f (x , u(x)) , Vx € D, 
\ u(x)\q£> = h(x)v (x), Vx 6 3D. 



E , 



(5.8) 



TD 



exp -2 / Jx(X°(s))ds+ q(X°(s))ds 



+ 



-(Vv)a(Vv)* -{b-b, Vv) ) (X°(s)) ds 



< CO, 



where E x indicates that the expectation is taken under P x . From Theo- 
rem 4.2, it follows that equation (5.7) admits a unique weak solution u. Set 
u{x) = h~ 1 (x)u(x). We will verify that u is a weak solution to equation (5.1). 

Indeed, for tp € W ' (D), since u(x) = h(x)u(x) is a weak solution to 
equation (5.7), it follows that 



»J=1 



2? l aij(x) 



d[h{x)u{x)\ d[h 1 (x)i/j] 



dxi 



dxj 



dx 



-E / Mx) - k(x) - (aVvUx^^^^h^ix^dx 
~[Jd OXi 

+ J (b — b, Vv (x))u(x)ip(x) dx 

(y v)a(V v)* (x)u(x)ij) dx — I q{x)u{x)ip{x) dx 



D 



f(x, u(x))ip(x) dx. 



Denote the terms on the left of the above equality, respectively, by T\, T-2, 
T 3 , T 4 , T 5 . Clearly, 

ri= if f a ..( x )^l^t d x-- V ( a-(x)^- — i>dx 
1_ 2 ^ D A) dxi dxj 2^ D A) dx, dx/ d 

i,j=l J ii3=± 

(5.9) +l£j^( x )^*t u ( x)das 

1,9 = 1 J 
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Using chain rules, rearranging terms, it turns out that 

T 2 + T 3 = - g JU*) ^ *** - g / d tog «(*) cfa 



d ■ 

(5.10) +V / + (oV«)* 

i=i ^ 



, d[Vm(x)] 
d 



/ (aVi))i(i)-r — ii(x)(ix + y^ / (aVu)j(a;) - — u(x)ipdx. 

"HJd oxi rrlJD oxi 



i=X " i=l 
In view of (5.4), 



(5.11) 



dx. 



Thus, 



T 2 + T 3 = - g ^ &,(*) 9 -^^dx-Y,J D k{x) g „(*) 
1 d . 

(5.12) +2E / ( aV ^( 3 
i=i ^ D 

1 (tin / ill) 

/ (aVu)i(i) - — u(:e) dec + / (aVv)Ax) —— u(x)ipdx. 

~{JD OXi 7^1 J D dx i 



d[i)u{x)] 
(x) — dx 

OXi 



After cancelations, it is now easy to see that 

d 



T 1 +T 2 + T 3 +T4 + T 5 = ^J2 J 

-E / w 

tl^ 



. du(x) dip , 

OXi OXj 



OXi 



(5.13) 

d f dtp f 
— 7 / b — — u(x)dx— / q(x)u(x)ip(x) dx 
;JD ox i JD 



f(x, u(x))i/j(x) dx. 



D 
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Since tp is arbitrary, we conclude that u is a weak solution of equation 
(5.1). Suppose u is a continuous weak solution to equation (5.1). Put u(x) = 
h{x)u(x). Reversing the above process, we see that u is a weak solution to 
equation (5.7). The uniqueness of the solution of equation (5.1) follows from 
that of equation (5.7). □ 
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